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For generic mesoscopic systems like quantum dots or nanoparticles, we study the Anderson orthogonality 
catastrophe (AOC) and Fermi edge singularities in photoabsorption spectra in a series of two papers. In the 
present paper we focus on AOC for a finite number of particles in discrete energy levels where, in contrast to 
the bulk situation, AOC is not complete. Moreover, fluctuations characteristic for mesoscopic systems lead to 
a broad distribution of AOC ground state overlaps. The fluctuations originate dominantly in the levels around 
the Fermi energy, and we derive an analytic expression for the probability distribution of AOC overlaps in the 
limit of strong perturbations. We address the formation of a bound state and its importance for symmetries 
between the overlap distributions for attractive and repulsive potentials. Our results are based on a random 
matrix model for the chaotic conduction electrons that are subject to a rank one perturbation corresponding, 
e.g., to the localized core hole generated in the photoabsorption process. 
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I. INTRODUCTION 

Among the simplest many-body problems one can consider 
is a Fermi sea of electrons interacting with a static impurity. 
Such an impurity is, for example, suddenly formed when a 
core or valence electron is excited above the Fermi energy by 
an incident photon in a bulk-metal photoabsorption or photo- 
luminescence experiment. Naively, one would expect a sharp 
onset of the absorption cross section at threshold with a signal 
proportional to the density of states in the conduction band 
above the Fermi energy. In experiments, deviations from this 
behavior were found in the form of peaked or rounded edges 
both in metals 1 and semiconductors. 2 These so-called Fermi 
edge singularities (FES) are a clear indication of the impor- 
tance of many-body processes that have been receiving en- 
during interest. 

For (clean) macroscopic systems, an essentially complete 
understanding of the mechanisms causing the Fermi edge sin- 
gularities has been obtained by the late eightiesi^i^ii In par- 
ticular, these works have emphasized the role of two com- 
peting effects. The first one is the Anderson orthogonality 
catastrophe (AOC), which is a property of the overlap be- 
tween many-body wavefunctions. However, in x-ray absorp- 
tion and photoabsorption processes, besides AOC there is an- 
other many-body response of the conduction electrons to the 
impurity potentialpiiiSi often referred to as Mahan's exciton, 
Mahan's enhancement, or the Mahan-Nozieres-DeDominicis 
(MND) contribution. It counteracts the effect of AOC in situa- 
tions where the dipole selections rules are fulfilled, and causes 
the photoabsorption cross section to diverge at the threshold. 

In the present series of two papers we are interested in these 
many -body effects in mesoscopic systemspiiiiii such as quan- 
tum dots or nanoparticles, and in particular in the changes 
with respect to the macroscopic (or bulk) case due to the con- 
fined geometry. In this first paper, we focus on the orthog- 
onality catastrophe mechanism for a finite system. We will 
more particularly consider the case of a contact (rank one), 



but not necessarily small, perturbation, with the localized per- 
turbation presented by the core hole left behind in photoab- 
sorption just being one example. In the second paper, that 
will be referred to as |lj in the following, we will discuss in 
detail Fermi edge singularities and mesoscopic photoabsorp- 
tion spectra that depend not only on AOC but crucially on the 
MND contribution. 

In a bulk-metal x-ray absorption experiment, the conduc- 
tion electrons respond to the abrupt, non-adiabatic appearance 
of a core hole generated upon photoexcitation of an inner shell 
electron into the conduction band by adjusting their single 
particle energies and wavefunctions slightly. Since the core 
hole represents an effectively positive charge the scattering 
potential will be attractive, and the energy levels will be low- 
ered. Although the overlap of the respective single particle 
states before and after the perturbation is very close to one, 
the many-body ground state overlap A tends to zero in the 
thermodynamic limit, A oc M~ E with M the number of con- 
duction electrons and the parameter e > 0. This effect, first 
noted by Anderson in 1967, 12 is the Anderson orthogonality 
catastrophe. As a result of AOC the cross section A(oj) for the 
absorption of photons with energy w will be power-law sup- 
pressed near threshold, explaining the above-mentioned Fermi 
edge singularity in the form of a rounded edge. 

AOC plays an important role in all processes where a fi- 
nite range scattering potential suddenly arises or changes. It 
is also prominent in the Mossbauer effect where AOC en- 
ables the recoilless emission or absorption of 7-rays by the 
whole crystal, but not by a single nucleus. 13 Another exam- 
ple is Kondo physics, except that in that case the impurity is 
dynamic, rather than static as considered here. 

Turning from macroscopic to mesoscopic systems opens a 
new venue to the study of AOC. To start with, the confined ge- 
ometry will naturally induce new behavior associated with the 
mesoscopic regime, among which are the lack of translational 
or rotational invariance, the existence of interference effects, 
and the related mesoscopic fluctuations. 91011 Furthermore the 
advances in the fabrication of mesoscopic and nanoscopic de- 
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vices provide a great flexibility in designing systems for which 
one can measure and control processes dominated by AOC. 
One example would correspond to the situation where an elec- 
tron tunnels through a system of two quantum dots joined via 
a small constriction. In that case the boundary conditions 
of the electrons on the far dot are changed once the tunnel- 
ing electron enters the first dot, 1415 and as a consequence a 
change of one electron appears as a sudden perturbation to the 
remaining electrons. Another example is a mesoscopic tun- 
nel junction containing a localized impurity level such that 
tunneling from the junction to the lead leaves behind a hole 
bound to the impurity. The resulting power-law singularities 
in the current- voltage characteristics were studied in Ref. llg 
and observed in experiments with two- and three-dimensional 
electrodes. 17 Large optical singularities in one-dimensional 
electron gases formed in semiconductor quantum wires were 
observed in Ref. [Hand theoretically explained as Fermi-edge 
singularities^ 

Aspects of AOC in disordered or mesoscopic systems have 
been addressed in various publicationsi22i2ii2£iS AOC in disor- 
dered simple metals and its influence on x-ray photoemission 
spectra (where the excited core electron leaves the metal and 
the edge behavior is determined by the AOC response alone) 
was studied by Chen and Kroha^ In Ref. |22j, Vallejos and 
coworkers studied AOC in chaotic mesoscopic systems for a 
weak perturbation modeled by a random matrix. They found 
large fluctuations of the AOC overlap that is non-zero in gen- 
eral due to the incomplete orthogonality. Small overlaps are 
correlated with avoided level crossings between the highest 
filled and the lowest empty level. This motivated the analytic 
study of a two-level model to which we refer below. Meso- 
scopic conductors with diffusive disorder were studied with 
respect to AOC by Gefen et alP The dependence of the av- 
erage AOC-overlap on the number M of electrons was found 
to be dimension dependent, and to show a stronger depen- 
dence on M in two dimensions than in higher dimensions. In 
agreement with the results by Vallejos et al. and the results 
we present here, a broad distribution of overlaps is seen, with 
the small-overlap tail of the distribution being related to few- 
level statistics near the Fermi-energy. Another key result was 
that the disorder-averaged AOC overlap depends non-mono- 
tonically on the disorder. This becomes clear in the limit of 
strong disorder where the electrons are strongly localized and 
cannot be sensitive to the addition of another impurity, and 
was confirmed in their numerical results. 

As already mentioned, here we will more specifically con- 
sider the case of a short range (rank one), not necessarily 
small, perturbationj^ii 2 ! which applies for the description of 
the perturbation created by the core hole. We will further- 
more consider a generic chaotic mesoscopic system, such as a 
quantum dot or nanoparticle, for which a random matrix the- 
ory (RMT) model of the conduction electrons applies. This 
particular regime will show some similarities with both the 
diffusive systems studied in Ref. 23 and the non-rank-one per- 
turbative regime investigated by Vallejos et al., 22 but will also 
display interesting specific properties. 

Beyond AOC, Fermi-edge singularities in the mesoscopic 
regime were considered in Refs.|8|,ll5|, andllo, We will con- 



sider the full mesoscopic x-ray edge problem and the resulting 
FES in photoabsorption spectra of quantum dots or nanopar- 
ticles in the second paper of this series, 11 and concentrate on 
AOC in the present, first paper. 

The paper is organized as follows. In Sec.[n]we will intro- 
duce our model of a rank one perturbation on discrete (bulk- 
like or chaotic) energy levels and apply it to the study of AOC. 
We investigate the properties of the ground state overlap prob- 
ability distribution P(|A| 2 ) and show that the phase shift at 
the Fermi energy determines the overlap fluctuations. The 
presence of a bound state is addressed in Sec. |lll| and its role 
for the symmetry of overlap distributions for positive and neg- 
ative perturbations is discussed. Knowing that the origin of 
the fluctuations lies in the levels nearest to the Fermi energy 
allows us to derive an analytic formula for the overlap distri- 
bution in the limit of strong perturbations in Sec. II VI We end 
this first paper with a summary on AOC in generic mesoscopic 
systems. 

II. ANDERSON ORTHOGONALITY CATASTROPHE FOR 
MESOSCOPIC CHAOTIC SYSTEMS 

Anderson orthogonality catastrophe refers to the vanishing 
of the overlap A between the system's many-body ground 
states before and after a finite range scattering potential is 
suddenly applied, and is complete only in the thermodynamic 
limitii^ In this Section we will discuss the effect of a finite 
number of particles in a mesoscopic system as well as the role 
of mesoscopic fluctuations. As the spin variable does not play 
a particular role here, we shall not take it into account (ex- 
cept for the discussion of the Friedel sum rule). We therefore 
describe the unperturbed system by the Hamiltonian 

JV-l 

Hp = e fe cj.c fc (1) 

k=0 

where the operator c k creates a particle in the unperturbed 
eigenstate yjfe(r) corresponding to the eigenenergy e&. We 
want to consider a situation where the perturbing potential is 
created by the hole left by a core electron excited into the con- 
duction band, which due to screening by the other electron is 
a short range potential. We shall therefore use for the pertur- 
bation a rank one, contact, potential 

V c = v c ilp{r c )f{r c ) , (2) 

[/(r c ) = ^ fc y>fc(r c )cfc], that only acts at the location r c of 
the core hole (fl is the volume in which the electrons are con- 
fined, and the parameter v c defines the strength of the poten- 
tial). Introducing c\ as creator of a particle in the perturbed 
orbital tp K (r) (we will use Greek indices to indicate the per- 
turbed system), we obtain the diagonal form of the perturbed 
Hamiltonian as 

H = H + V c = X ^c K ■ (3) 

K 

We will refer to the many-body, Slater-determinant ground 
states of the unperturbed and the perturbed system as |$o) 
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and l^o); respectively. Their overlap is A = (^ol^o)- In the 
ground state the M conduction electrons fill the lowest of the 
N levels, i.e. the levels i — ... M — 1 [see also Fig. 0a)]. 
Whenever practical, we will use the index i (/i for perturbed 
levels) for levels below Ep and j (7) for levels above Ep, and 
k (k) when reference to all levels is made. 



A. Bulk-like systems 

In the bulk the orbital momentum I is a good quantum num- 
ber, and, because the core hole potential (|2j is rotationally 
invariant, only the I = sector couples to the perturbation. 
Within this sector, the unperturbed energy levels {ek} near 
the Fermi energy can be taken equidistantly spaced and, with 
a proper choice of the origin of the phase of the wavefunctions 
we can take </>fc(r c ) = const. = 1/vfi. 

In the following, we will call "bulklike" a model which, as 
the I = sector of a bulk system, has no fluctuations for the 
eigenlevel spacing nor for the amplitude of the wavefunction 
at the impurity. We shall, however, assume finite the number 
of levels within the bandwidth. In this subsection we shall 
review briefly some properties of this bulklike system. 

In the absence of mesoscopic fluctuations, the perturbation 

is 
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Introducing the transformation matrix a = (cifc K ), 



N-l 

E 

fc=0 



(4) 



one can check easily that the perturbed eigenstates {X K } and 



the coefficients akn fulfill the familiar equations 
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= l-v c 



N-l 
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(5) 
(6) 

(7) 



Eq. (|5} allows one to find the perturbed eigenvalues A K as a 
function of the unperturbed energies and the perturbation 
strength v c . This is illustrated in Fig.[^where the perturbed 
eigenvalues A K are found as intersection points of the function 



A - e k 



with the horizontal line l/v c . We have assumed a negative 
perturbation v c < to account for the attractive perturbation 
caused by the core hole potential. 

The effect of the perturbation and the shift between per- 
turbed and unperturbed levels is generally characterized by 
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FIG. 1: (color online) Perturbation of a bulklike material with six 
initially equidistant energy levels (0 < < 2n, mean level spac- 
ing d) by a rank one perturbation of strength v c < 0. The perturbed 
eigenvalues A re arise from a graphical solution of Eq. J5} and are 
shifted to lower energies (attractive interaction). This shift increases 
as the perturbation becomes stronger; see the horizontal lines corre- 
sponding to two different potential strengths. The strongest perturba- 
tion shifts the perturbed levels right to the middle of two neighboring 
unperturbed levels e^, 6k+i (strictly, this applies in the center of the 
band - see the discussion of boundary effects). 



the partial-wave phase shifts 61 at the Fermi energy for the 
orbital channel I. The phase shifts obey the Friedel sum rule 



Z = ^2(2Z + l)V7r 



(8) 



where the factor of two accounts for spin and Z = — 1 in the 
case of photoabsorption. For a localized rank-one perturba- 
tion which, as mentioned above, is spherically symmetric, all 
61 with I are zero. In this case the Friedel sum rule Eq. (|8j 
implies 60 = —tt/2. 

When Ep is in the middle of the conduction band, So is 
related to the potential v c and the mean level spacing d by 21 



Sq = arctan 



irv c 



(9) 



Sq = —tt/2 therefore corresponds to a strong perturbation 
\v c /d\^>l. We will therefore be mainly interested in AOC in 
this limit, although we shall also consider other cases. 

For the finite number of particles that we consider in the fol- 
lowing, boundary effects are important. First, the phase shift 
varies (monotonically) with the level number, and only for 
levels at the band center Eq. l|9} holds. For levels off-center, 
the number of energy levels on the left and right is unbalanced 
and leads to additional pushing from the side with the surfeit. 
This can already be guessed from inspecting Fig.0 

This effect can easily be taken into account by noting that, 
for any perturbed level X K [assuming for instance < k < 
(N — l)/2], Eq. l|5} can we rewritten as 

2k- 1 

y =- 

f-^ X K - E k V K 

K = (j 

where we have defined 
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Since k is now in the middle of the range [0, 2k — 1], we see 
that Eq. can be generalized to give a level-dependent phase 
shift 6 K = arctan7ri; K /d, where evaluating the sum as an in- 
tegral yields 
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if k e o, 



N — 1 



(12) 
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N - 1 



,N- 1 



Note that away from the band center, 6 K is no longer confined 
to the interval [— 7r/2,7r/2] but can take values in the range 

[-7T,7r]. 

A second boundary effect becomes important for strong 
perturbations v c j d < — 1 that induce a large shift of the low- 
est perturbed energy level towards — oo (see Fig.Q. We will 
discuss the effect of such a bound state in detail in Sec.lIIII 

A remarkable property of rank one perturbations such as a 
contact potential is that the overlap between the many-body 
ground states |$ ) an d l*o) of H and H can be expressed 



M-l N- 



|A| 2 = |<*o|$o) 



ji n ft ~ ?w c . ~ Xi l ■ " 3 » 



11 (A, 



In other words, all the required information is contained in 
the eigenvalues {e}, {A}. In the bulklike situation we have 
considered so far, the {A} are uniquely determined from the 
{e}, yielding the bulklike overlap |A;,| 2 . Analytic formulas 
for |A b | 2 were discussed in Ref.0 For half filling, M = N/2, 
and assuming a level-independent phase shift, the relation 



[ '-r 



(14) 



is easily derived. The factor T ~ | SfcLi w ~ 0.822 is a 
correction to the exponent due to the discreteness of the levels. 



B. Mesoscopic systems 

We now consider a generic chaotic mesoscopic system. In 
that case, because of the confinement, both energies and wave- 
functions will display mesoscopic fluctuations from level to 
level or as an external parameter is variedi^iiSiii This implies 
that the effective strength of the contact potential felt by elec- 
tron k cannot be assumed constant. More precisely, introduc- 
ing Ufe = \/Titpk(r c ), such that (|iifc| 2 ) = 1, we now have 

t. 



V C =V C 2_ J u* k u k 'c\ck 



kk' 



Accordingly, the generalized version of Eqs. (|5}-(Q f° r me 
perturbed energy levels \ K reads 
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k=0 



(15) 
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tN 



VCR 

^ dknUk 
k=0 



(16) 
(17) 

(18) 

(19) 



and the bulklike situation is recovered by setting u 2 . = 
(|ufe| 2 ) = 1. Note that we can choose the phase in the wave- 
functions ipk and ipK such that their value at r c is real. It is 
then evident that u k , a kK , and u K are also real. 

For a chaotic system, a good model for the eigenvalues 
and eigenvector fluctuations is given by the classic ensembles 
of Random Matrix Theory (RMT). 24 ' 25 Most commonly the 
Gaussian orthogonal (unitary) ensembles, GOE (GUE), are 
used to describe the energy eigenvalues of chaotic systems in 
the presence (absence) of time reversal symmetry. Here how- 
ever, because we would like to model systems with constant 
mean density of states, we will rather consider the correspond- 
ing Circular Ensembles COE and CUE 25 which display the 
same fluctuations but have by construction a mean density of 
eigenphases exp(iefe) that is uniformly distributed on the unit 
circle. Since time reversal symmetry can easily be broken by 
applying a magnetic field, we will address both ensembles. 

The fluctuations in the {e} and {A} are not independent. 
Under the hypothesis that the unperturbed eigenenergies are 
distributed according to the circular ensembles, and accord- 
ingly that the eigenvectors follow a Porter- Thomas distribu- 
tion, the joint probability distribution P({e}, {A}) was shown 
by Aleiner and Matveev 21 to be 



P({e},{A})oc 



n 



1,3 



A/ 



-0/2 



Aj) _£ 

e 2 



£ 4 (Ai-ei)M 



(20) 

with the constraint e,_i < Aj < ej and j3 = 1 in the COE case 
whereas (3 = 2 for CUE. In the middle of the band Vi = v c , 
away from the band center, boundary effects have to be taken 
into account according to Eq. ill 2b . It is interesting that a joint 
probability distribution can also be found in the more general, 
non-rank-one case 26 that we do not consider here. 

The ground state overlap A is expressed entirely in terms 
of the unperturbed and perturbed eigenvalues {e} and {A} 
[see Eq. ( I13H . Therefore, the joint distribution i20\ contains 
all the information required to build the overlap distribution 
P(|A| 2 ). 

In Section IIVI we shall see that a slightly different route 
is actually more effective to obtain analytic expressions for 
P(\ A| 2 ). However, the present approach is perfectly adapted 
to constructing numerically this distribution for any value of 
the parameters M, N, v c , and 0. Indeed, we can use Eq. J20l > 
as the basis of a simple Metropolis algorithm to generate en- 
sembles of ({e}, {A}) with the proper joint probability distri- 
bution. Applying Eq. i\3l . we then build the distribution of 
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FIG. 2: (color online) Distribution of ground-state overlaps in a 
mesoscopic system as the number of particles is increased at half- 
filling, M = N/2, for the (a) CUE and (b) COE cases. The pertur- 
bation strength is chosen intermediate, v c /d — —0.5. The dash- 
dotted lines indicate the bulklike overlaps |A(,| 2 . The curves il- 
lustrate the power-law characteristic of AOC that makes it a truly 
thermodynamic-limit effect: even with M = 1000 electrons in the 
system, the most probable overlap |A| 2 is as big as 0.46, indicat- 
ing that a bulk-based discussion of AOC in mesoscopic systems can 
be inaccurate. Note that for this intermediate perturbation strength 
the overlap distribution exhibits a clear maximum and a difference 
in CUE vs. COE for small overlaps. (In order to make the statistical 
error essentially negligible, 20000 realizations of the random matrix 
were used in producing each curve.) 



ground state overlaps. Examples for overlap distributions ob- 
tained in this way are shown in Fig. [2] for systems with and 
without time-reversal symmetry and an intermediate pertur- 
bation strength. The evolution of AOC as the number of par- 
ticles is increased is clearly visible as a shift of the probabil- 
ity density to smaller values of |A| 2 . However, since AOC 
is a power-law effect even TV = 1000 is far away from the 
thermodynamic limit, cf. Eq.( ll4t that can also be used as a 
rough estimate for the mean overlap in the mesoscopic case. 
We point out the different behavior for very small overlaps in 
the CUE vs. COE case. It originates in the difference of the 
Porter- Thomas distributions and will be discussed in Sec. llVI 

In the metallic x-ray edge problem all relevant properties 
are known to depend only on the phase shift at the Fermi en- 
ergy, and a natural question is whether this remains true in 
the mesoscopic case for the groundstate overlap. In Fig.|5Ja) 
we therefore compare overlap distributions for different pa- 
rameter sets {N, M,v c /d} all of which yield a phase shift 
Sf ~ Bictaxi(TTVi—M / d) ~ — 7r/2 at the Fermi energy Ep. 
The curves differ, as do the corresponding bulklike overlaps 
|Af,| 2 indicated by the vertical lines. However, after scaling 
by the bulklike overlaps all curves coincide, Fig.|3jb). There- 
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FIG. 3: (color online) Ground state overlap in the absence of time- 
reversal symmetry (CUE case) for different parameter realizations 
corresponding to a phase shift of approximately —n/2 at the Fermi 
energy, (a) Before, (b) after scaling with the bulklike result |A;,| 2 
(dash-dotted lines, |Af,| 2 increases with the legend entries). The 
scaled result for COE is given in the inset of (b), see Ref. 8 for the 
unsealed result in that case. 

fore, whereas the individual P(|A| 2 ) depend independently 
on N, M, and V c /d, the fluctuations of the overlap depend on 
the value of Sf alone. 

This fact allows one to provide the overlap distribution as a 
universal curve for a given phase shift Sf- In Figs.|3a) and 
|5ja) results are shown for N = 100 levels subject to a strong 
perturbation. The phase shift Sf is varied by increasing the 
filling of the band from 2 to 98 electrons. The distributions are 
not scaled by the bulk overlap in order to allow easier compar- 
ison of the curves when they all range from to 1 . 



C. Fluctuations of levels beyond the Thouless energy 

Strictly speaking a random matrix model applies only to 
chaotic energy levels on a scale of the Thouless energy 
2J Th 4£iiii2iS In the semiclassical limit £Jrh is much larger 
than the mean level spacing, but significantly smaller than the 
Fermi energy Ep- Therefore one may question the model we 
use where the full bandwidth is described by RMT. It turns 
out, however, that the overlap distributions depend only on 
the fluctuation of the levels near the Fermi energy. In order 
to demonstrate this, we introduce a range-n approximation in 
which only n of the levels on either side of the Fermi energy 
are left free to fluctuate, while all the others are kept fixed at 
their mean position. To gauge the accuracy of this approxima- 
tion, we perform a Kolmogorov-Smirnov test, 2 - In Fig.[6]the 
result of this test, namely the maximum deviation of the range- 
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FIG. 4: (color online) (a) Probability distribution P(|A| 2 ) for the 
ground state overlap (COE, N = 100,v c /d = —10). The number 
M of electrons is varied, corresponding to a varying phase shift Sp 
at the Fermi energy. The curves (here not scaled by the bulklike over- 
laps) illustrate the characteristics of P(| A| 2 ) for a given phase shift. 
For comparison, recall that at constant filling, a decreasing phase 
shift corresponds to increasing perturbation strength, (b) Compari- 
son of the bulk-scaled probability distributions of the ground state 
overlaps for a negative (solid line, |A| 2 ) and positive (dashed line) 
perturbation. The overlap |A + | 2 is, up to the negligible influence of 
the run-away states, equivalent to the overlap | Aom | 2 shown here as 
dashed line. It will be of importance in the discussions of the x-ray 
absorption part. The pairing of curves whose phase shifts moduli add 
up to 7r is nicely illustrated in this graph. 



FIG. 5: (color online) Same as in Fig. |4| but for CUE. Note the 
difference in the distributions for small overlaps | A| . 



These so-called run-away levels correspond to the formation 
of a bound statei^LS which entirely screens the impurity. 



They can be found, respectively, at 



Ao — eo — — 
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n cumulative distribution from the exact cumulative distribu- 
tion of the overlap, is shown. Clearly, the deviation decreases 
to almost zero as soon as the fluctuations of a few levels near 
Ep are included. The convergence with n appears further- 
more faster for the larger perturbation strengths in which we 
are interested for the x-ray edge problem. This retrospectively 
justifies the use of a full random matrix to model the statisti- 
cal properties of our chaotic systems, as the levels beyond the 
Thouless energy for which such a description does not hold 
do not affect the overlap distribution. 



III. BOUND STATES AND SYMMETRY BETWEEN 
POSITIVE VS. NEGATIVE PERTURBATION 

As soon as the perturbation strength becomes significant, 
i.e. for attractive perturbations v c /d < — 1, the shift of 
the lowest perturbed level, Ao — en, starts to become large. 
Analogously, for positive perturbation strengths larger than 
v c /d > 1, the same boundary effect now shifts the highest 
perturbed level A^v-i to very high energies, see Figs.^ an d0 



The fluctuations of these levels will be neglected, since in 
the limit of very strong perturbations, |u c |/d 3> 1, the run- 
away level does not affect the value for the overlap, Eq. dl 3I >. 
because the terms involving Ao (or Ajy-i, respectively) yield 
ratios of one. Therefore, one has only to consider the N 
unperturbed levels {e} and the N — 1 perturbed levels A K 
"in between". In this case of very strong perturbations, the 
perturbed levels are maximally shifted: in the center of the 
band, the X K lie, on average, exactly in the middle between 
the two neighboring unperturbed levels e K _i and e K (or e K 
and e K+ i). Therefore, a strong negative perturbation v c (i.e. 
v c /d -C — 1, Sp = — 7r/2) and a strong positive perturba- 
tion u+ (i.e. v£/d 3> 1, Sp = 7r/2) are indistinguishable, 
see Fig. This means in particular that the corresponding 
overlap distributions P(|A| 2 ) and P(|A + | 2 ) are identical. 

In the case of intermediate perturbation strengths v c /d ~ 
±1, a level shift Sp associated with the negative perturbation 
corresponds to a phase shift Sp = it — \Sp\ when the level 
scheme is interpreted as originating from some positive per- 
turbation u+ ( 7^ |w c | in general), see Fig.Qc). We therefore 
expect again the two corresponding bulk-scaled overlap dis- 
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FIG. 6: Accuracy of the range-n approximation using the 
Kolmogorov-Smirnov test (COE, N = 100, M = 50). The maxi- 
mum deviation of the range-n cumulative distribution from the exact 
cumulative distribution of the overlap is shown as a function of n. 
The rapid approach to 0, especially in the strong perturbation case, 
shows that fluctuation of only those levels near Ef is important. 



tributions to be identical, and the (symmetry) relation 

|A| 2 \ „/|A+' 2 ' 

v c ,8 f 



P 



|A^I 2 



(21) 



-IM 



to hold. 

In Figs. |4] and [5] overlap distributions for phase shifts > 
Sf > — 7T are shown. Starting with an almost full band 
{M = 98 particles on N = 100 levels) and a small nega- 
tive phase shift, Sp is lowered (the effective potential strength 
raised) as the filling is decreased, passing Sf ~ — 7r/2 at half 
filling. The corresponding overlap distributions are shown in 
Fig-Ek for the COE situation and in Fig. [5^ for a CUE case. 
In the lower panels (b) of Figs. |4] and [5] overlap distributions 
belonging to a negative phase shift Sf (and v c < 0) and the as- 



sociated positive phase shift Si 



\Sf\ (and > 0) are 



compared and seen to (almost) coincide pairwise after scaling 
by the respective bulk overlaps (the small deviations are due 
to the finite number of discrete levels available). Note that the 
bulk overlaps |AjJ~| 2 are very similar in each pair, but rather 
different for different pairs, thereby causing the stretching of 
the x-axis. 

In the context of the x-ray edge problem (strong negative 
perturbation), a situation that is equivalent to a positive per- 
turbation becomes of importance when the core electron is 
excited to the level Am just above the Fermi energy Ep ■ The 
consequences and more details will be discussed in|lj 



IV. ANALYTIC RESULTS FOR THE OVERLAP 
DISTRIBUTION 

We finish with the derivation of analytic results for the over- 
lap distribution. We argued in Sec. lIICI that the fluctuations in 
the ground state overlap are dominated by fluctuations of the 
energy levels around the Fermi energy (see again Fig.|6j. We 
now want to deepen this discussion and investigate to what 
extent an analytic understanding of overlap probability dis- 
tributions P(|A| 2 ) is possible. In particular we want to un- 



N-1- 



(a) i/ c «-1 



J d — e- 



(b) V c »1 



(C) |VJ ~ 1 



M-r 




FIG. 7: (color online) Run-away level for strong perturbation which 
is (a) positive or (b) negative. For simplicity, the bulklike situation 
is shown. In the case of intermediate perturbations (c), the perturbed 
levels \ K are not shifted to the center of the unperturbed level inter- 
val efc . The phase shifts Sp (negative perturbation) and 5 J (positive 
perturbation) are related by \8^\ + 8p~ — ir. In the case of strong 
negative perturbations, as in the x-ray problem, Ao forms a bound 
state that screens the perturbation. 



derstand two characteristic qualitative features of the overlap 
distribution (see for example Fig.[3j: firstly P(|A| 2 = 0) is 
finite and characteristically different for the COE vs. CUE sit- 
uation, and secondly P(|A| 2 = 1) is zero, i.e. P(|A| 2 ) drops 
to zero before the natural bound | A| 2 = 1 is reached. 

To begin, remember that in Eq. dl 31 for the overlap |A| 2 , 
all energy differences are taken between empty (j > M) and 
occupied (i < M — 1) levels. The smallest differences (of the 
order of the mean level spacing d) contain the energy levels 
£a/, £m-i> and Aa/- We expect therefore that the small A 
part of the overlap distribution will be mainly given by the 
fluctuation of these terms. 

To gain an analytic understanding of the overlap distribu- 
tion, we therefore consider only the level closest to the Fermi 
energy as fluctuating whereas all other levels are assumed at 
their bulklike values (range- 1 approximation). We therefore 
want to evaluate fluctuations of (Am — £m-i)/(em — £m-i)- 
For simplicity, we consider the half-filled case. Note that a 
similar two-level problem was studied by Vallejos et al. in 
Ref. |22j, motivated by their observation that very small AOC 
overlaps were related to avoided level-crossings at the Fermi 
energy. In our present model, Eq. dl 51 for the rank one per- 
turbation determines the position of the perturbed levels as 
solutions of the equation 



A - e k 



d 



(22) 



where A k = \u k \ 
function intensities, 

P(A) = 



for CUE 

e -A/2 for C q E 



(23) 



[{A) = 1 and (A) 2 = 2(3) for CUE (COE)]. We will take 
only the fluctuations of the amplitudes Am-i and Am for 
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the wavefunctions <\>m-\ and 4>u explicitly into account. Ac- 
cordingly, we separate the first two terms of the sum and com- 
prise the remaining terms in a function Sf that is assumed to 
be a Gaussian random variable and contains the joint effect of 
all other levels and intensities, 



/(A) = 



dA 



M-l 



dA 



M 



A — cm-i A — ejvf 



(24) 



Equation 1241 makes it possible already to understand at 
a qualitative level the finite probability for having |A| 2 = 
which corresponds to having Am — ejwr-i = 0. How- 
ever, setting A = Am in Eq. J24l > we see that the ratio 
Am-i/(Am — £m-i) has to be finite. We conclude that hav- 
ing Am — £m-i = is related to having Am-i = 0. From 
Eq. O we find P(A = 0) = 1 [or oo] in the CUE [COE] 
case, thereby explaining both P(\ A| 2 = 0) > and the peak 
at |A| 2 = 0intheCOEcase. 

To obtain a more quantitative description, we investigate 
Sf more closely. According to the assumptions of the range- 
1 approximation, we replace the energy levels in Sf by their 
bulklike values. We furthermore neglect boundary effects that 
would make the phase shift dependent on the level number 
which is a good approximation for N — > oo at half-filling. 
Fixing the energy origin in the middle of the mean position of 
e&f_i and em, we express Sf in the symmetric form 



.4 



k+M 



A_ 



(fe+l)+M 



fc=l 



(25) 



Fluctuations in Sf are due to the fluctuations in the intensi- 
ties A±k that follow the Porter- Thomas distributions given in 
Eq. (I23i . Expanding the mean value (Sf) about the center of 
the band, A = 0, yields in lowest order (Sf) — ~2XSi with 



Si = 



E 



i 



1)2 



0.9348 . 



(26) 



Comprising all the randomness of the amplitudes A)-^m-1,M 
in one random Gaussian variable^ with P(£) = / 2 /V2n, 
and writing explicitly the scale \aiSf of the fluctuations 
around the mean value, we arrive at the approximation 



Sf = (Sf) + vW/ £ = -2XS, + Vvar^/ £ (27) 



with 



var<5/ = 



2Si for CUE 
4Si for COE. 



(28) 



Setting £m — £j/-i = s, the probability distribution P(A|s) 
is obtained by integrating over the distributions of Am- i , Am 
and £. These integrals can only partially be solved analytically 
and, leaving out the further details of the calculation, the result 

for the CUE case (p — au — aM-i,&o — Am-i/(A + 
s/2),ai = -A M /(X-s/2))is 



P(X\s) 



1 



2s\JitSi 
and for the COE case 
P(\\s) = 



d/i < exp 



/LtA 



™2. 



exp 



dfx 



exp 



8Si V 



25iA- 



(29) 
(30) 



2 

25i - A 2 ) + /iA 



T 



£/i 

4 



with A'o, K\ the modified Bessel functions. 

The remaining integral over /j, can be performed in the CUE case using the (inverse) error function. Splitting the //-integration 
into integration from to oo and from — oo to 0, yielding results P + (\\s) and P_(A|s), and using the substitutions A = 

2 Si A + l/v c , B = 25i (s 2 /4- A 2 ), a= l/(45i), and fe± = [s ± l/(S*iu c )]/2 we find 

exp(-^4 2 a) 



P±(A| S ) = 



2aC 



± 



exp 



erfc 



2^/E 



(31) 



and their sum yields P(A|s) = P+(A|s) + P_(A|s). 

This expression simplifies considerably in the limit of very strong perturbations, \v c \ ^> 1. Then, the result becomes indepen- 
dent of v c (in particular independent of the sign - very strong attractive and repulsive interactions are equivalent), and takes the 
form 



Poo(A|s) 



exp(-A 2 Si) 



sVSi + Vtt (1 - 2A 2 Si) exp 
I 



S 2 5! 



erfc 



(32) 



Clearly, this probability distribution is symmetric with respect to A = as it should be. 
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The contribution to the overlap, Eq. d 1 3 b . that contains the 
energy levels closest to the Fermi energy is represented by 
(Am - e M -i)/(eM - e M -x) = (A + s/2)/s in our present 
notation. To obtain the probability distribution P F 1 (A) of this 
quantity, we have to take into account the distribution of near- 
est neighbor spacings s = eu—^M-i according to the Wigner 
surmise, 



:S2 



P(S) = 



s 2 exp(-|s 2 ) for CUE 
for COE . 



f s exp (-|{ 



(33) 



The final result then reads 



POO 

P F1 {\)= ds sP(X\s)P(s) , (34) 
Jo 

This integrated result is however essentially indistinguishable 
fromP(A|s = 1). 

The probability distribution of the range- 1 approxima- 
tion, P rl (|A rl | 2 ), is straightforwardly obtained from P F1 (\) 
by scaling with the bulklike quantity |A;,| 2 /(1 — Sf/it) 
[cf. Eq. ( I35> below] and shifting the argument by 1/2. The 
resulting curves are presented in Fig.[8]for weak, intermedi- 
ate, and strong perturbations. This analytic range- 1 result is 
compared to the full overlap distribution P(|A| 2 ) as well as 
the exact range- 1 and range-2 approximations 34 (CUE case, 
see Ref. 8 for the corresponding COE figure). The agreement 
between the analytic approximation and the exact range- 1 re- 
sults is excellent. The analytic range- 1 result provides in par- 
ticular a very good estimate for the full overlap distribution 
at small overlaps, independent of the perturbation strength. 
Moreover, in the limit of strong perturbations in which we 
are especially interested for the x-ray edge problem, it pro- 
vides a reasonable estimate of the full probability distribution 
P( | A | 2 ) . The increase in the degree of symmetry of the range- 
1 result (with respect to the maximum of the curve) as the per- 
turbation strength is increased is clearly visible and confirmed 
by Eq. OH- 

We would like to mention that Vallejos and coworkers in 
Ref. I22L and, based on their work, Gefen and coworkers in 
Ref.|22l found similar results, in particular that the distribution 
of small-overlaps is well represented by a two-level analytic 
model. However, there are few differences. For example, the 
difference between the (Gaussian) orthogonal and unitary case 
was less pronounced in their models, especially the increase of 
the probability P(|A| 2 ) towards |A| 2 — > 0, characteristic of 
our model, was not seen. Further, simple analytic expressions 
similar to the expressions Eqs. J29M 30) could only be found 
for the orthogonal situation in Ref. |22j 
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FIG. 8: (color online) Probability distributions P(|A| 2 ) for the 
ground state overlap (N — 100, M = 50, CUE) for increas- 
ing potential strength: (a) weak (v c /d = —0.25), (b) intermediate 
{v c /d = —1), and (c) strong (v c /d — —10) perturbation. Shown 
are the exact distribution (dots) as well as the range- 1 (triangles) and 
range-2 (diamonds) approximations. Clearly, the fluctuations come 
from the levels closest to the Fermi energy, and the range- 1 approx- 
imation provides a good estimate already for intermediate perturba- 
tion strengths as in (b). The crosses denote the analytic result for the 
range- 1 approximation, showing excellent agreement with the nu- 
merical results. 



smaller than 1 and, furthermore, get smaller as the perturba- 
tion gets stronger. Accordingly, the upper boundary 1 of the 



first term is shifted to a value |A' 



< 1 that decreases 



with increasing perturbation strength as is clearly visible in 
Fig-El The value | A^ ax | 2 provides a good estimate for the up- 
per boundary of the overlap | A| 2 , although there are of course 
corrections from higher order approximations that smear the 



Eventually, we briefly remark on having P(|A| 2 = 1) = (total) over l a p distribution P(|A| 2 ) about |A^ ax | 2 , see for ex. 



for finite perturbation strength. Within a range- 1 -like approx 
imation scheme, Eq. Jl 31 can be written as 



rl|2 _ (Am - £M-l) d_ 

1 ~ (Am-Am-i) d (!-£)' 



IAJ 2 



(35) 



where we have normalized all differences by their bulklike 
values. Only the first term is fluctuating, between and 
1, whereas the last two terms are numbers that are strictly 



ample Fig.[8lb). 



V. SUMMARY AND CONCLUSIONS 

In this paper we have studied the Anderson orthogonality 
catastrophe in generic mesoscopic systems with a finite num- 
ber of electrons in discrete energy levels. Motivated by the 
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study of mesoscopic photoabsorption spectra that will be the 
subject of the second paper II in this series, we have used the 
model of a rank one perturbation describing the perturbation 
left behind after the photoexcitation of a core electron into the 
conduction band. However, the basic underlying physics is 
much more general, and the model can be applied to any sit- 
uation where a sudden, localized (rank one) perturbation acts 
on a finite number of chaotic electrons. 

As a result of such a sudden perturbation, the overlap 
between the unperturbed and perturbed many-body ground 
states vanishes in the thermodynamic limit (AOC). In meso- 
scopic systems, even with few thousand electrons, this effect 
will be incomplete. Fluctuations cause a broad distribution 
of overlaps that is, however, bounded from above by a value 
smaller than 1 (for not too small perturbations) whereas there 
is a non-zero probability for having zero-overlap. The distri- 
bution of small overlaps is distinctly different for systems with 
or without time reversal symmetry (COE vs. CUE) and traces 



back to the difference in the Porter- Thomas distribution of 
their single-particle wavefunction intensities. The fluctuations 
originate from the levels around the Fermi energy (within a 
range smaller than the Thouless energy for systems we are 
interested in), and we derived an analytic expression for the 
lowest order (range 1) result that reproduces the distribution 
of small overlaps and is a good approximation for the whole 
distribution in case of strong perturbations. 
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